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Abstract 

We present the next-to-leading order results for universal non-forward anomalous dimensions 
of Wilson twist- 2 operators in N=4 super symmetric Yang-Mills theory. The whole calculation 
was performed using supersymmetric Ward identities derived in this paper together with already 
known QCD results and does not involve any additional calculation of diagrams. We also 
considered one particular limit of our result, which could potentially be interesting in the context 
O of AdS/CFT correspondence. 

^ 1 Introduction 

| Parton distributions in QCD satisfy the Balitsky-Fadin-Kuraev-Lipatov (BFKL) pQ and Dokshitzer- 
5-h ; Gribov-Lipatov-Altarelli-Parizi (DGLAP) [2] equations. Up to day even next-to-leading corrections 
to these equations are known jSJIl]. Moreover, recently, similar results were also obtained for su- 
persymmetric theories El ■ An idea to consider these equations in the case of supersymmetric 
theories is based on common expectations, that the presence of extra high symmetry may significantly 
simplify them as well as their analysis. For example, it was already known, that all quasi-partonic 
operators in N = 1 SYM form supermultiplet of operators jSJ, having the same single universal 
anomalous dimension with shifted argument, which at that time was computed in leading order 
(LO) of perturbation theory. Calculations in N = 4 SYM gave even more remarkable results - the 
eigenvalues of the integral kernels in the evolution equations for quasi-partonic operators are propor- 
tional to — 1) — ^(l) [HJE]], which means that these evolution equations in the multicolour limit 
are equivalent to the Schrodinger equation for the integrable Heisenberg spin model similar to the 
one found in the Regge limit Moreover, it was shown, that in the maximally supersymmetric 
N = 4 Yang-Mills theory there is a deep relation between BFKL and DGLAP evolution equations 
In particular, the anomalous dimensions of Wilson twist-2 operators in N = 4 SYM could be found 
from the eigenvalues of the kernel of BFKL equation. The corresponding next-to-leading order (NLO) 
calculations showed, that many of these findings are valid also in higher orders of perturbation theory. 
Moreover, as in leading order they most fully realize in maximally supersymmetric N = 4 SYM. 
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However, some of the properties of these equations, valid at leading order, are violated at higher 
orders. For example, the conformal invariance of the theory in leading order allows us to construct 
multiplicatively renormalized quasi-partonic operators |H] up to this order in perturbation theory. 
But, in next-to- leading order in perturbation theory multiplicative renormalization of conformal 
operators is violated due to necessity in regularization of arising ultraviolet divergences, which is 
responsible for violation of conformal symmetry. 

The Efremov-Radyushkin-Brodsky-Lepage (ER-BL) equation jT5j may be viewed as some kind 
of the generalization of the DGLAP equation for the case of non-forward distribution functions. 
In this case in and out hadronic states in matrix elements of operators are different, what allows 
us to study scale properties of hadron wave functions. In the latter case the matrix elements of 
corresponding operators are considered between vacuum and hadronic state. In leading order the 
kernel of ER-BL evolution equation, after Mellin transformation was performed, is diagonal in basis 
of Gegenbauer polynomials [T51 ll4|. This is equivalent to the fact, that multiplicatively renormalized 
operators of ER-BL evolution kernel coincide in leading order with conformal operators constructed 
from Gegenbauer polynomials |15j . The direct calculation of next-to-leading order corrections to 
the evolution kernel of this equation turned out to be much more difficult problem compared to the 
forward case. It was first performed for the case of non-singlet quark operator in EI], where 
in the last reference an advanced computation method was developed and some analysis of general 
properties of evolution kernel and corresponding anomalous dimension matrix was given. In this 
case as well as in the case of model a solution of ER-BL evolution equation was constructed 
at next-to-leading order of perturbation theory [T7]. It was noticed, that at next-to-leading order 
conformal invariance is violated and anomalous dimension matrix develops non-diagonal part in the 
basis of conformal operators. 

The source of conformal symmetry breaking was identified later in |T%1 119j. It was found, that 
non-diagonal part of anomalous dimension matrix arises entirely due to the violation of special con- 
formal symmetry. Moreover, a framework based on the analysis of broken conformal Ward identities 
was proposed, which allowed to obtain for the first time next-to-leading order corrections to the 
nonforward anomalous dimension matrix in singlet case both in QED [20] and in QCD [21]. The 
generalization of these results for the case of supersymmetric theories involves a number of relations 
between elements of anomalous dimensions matrix [22], some of which were already known from the 
analysis of forward limit (last paper in Ref. [2],jH])- Later, the use of these supersymmetry relations, 
which could be also easily written for the corresponding evolution kernel matrix [HI [22], allowed to 
determine for the first time the QCD ER-BL evolution kernels in the singlet case both for odd [23] 
and even [21] parity distribution amplitudes (see also j2H])- 

In this paper we, using already known results from QCD, found closed analytical expression 
for universal non-forward anomalous dimension for maximally supersymmetric N = 4 Yang-Mills 
theory. Earlier, the authors of [H] used already known QCD results jl] to obtain expression for NLO 
universal forward anomalous dimension of Wilson twist-2 operators in N = 4 SYM, which were later 
confirmed by direct calculations in Ref. [7j. A problem, we are solving here, may be formulated as 
follows: how knowing QCD results one may obtain analogous results in supersymmetric theory with 
minimum efforts. 

It turns out, that owing to remarkable properties of N = 4 SYM it is possible to derive an 
expression for universal non-forward anomalous dimension without any additional calculation beyond 
those already done for QCD. 

Our result could be also interesting in the context of AdS/CFT correspondence [20] ■ Namely, 
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there are some calculations of the anomalous dimensions of such operators in the limit of large j 
(Lorentz spin) from both sides of the AdS/CFT correspondence j23 123 EH EDI HI EI]- There are, 
also some predictions for anomalous dimensions of other types of operators corresponding to multi- 
spin solutions in AdS^ x 5*5 space from string theory side [32] \ partially confirmed by field theory 
calculations. 

It should be noted, that up to this moment only diagonal part of anomalous dimension matrix 
have been studied in the context of AdS/CFT correspondence and it would be interesting to compare 
the results for its non-diagonal part with the appropriate result from string theory. 

The present paper is organized as follows. First in next two subsections we define the lagrangian 
of TV = 4 SYM and introduce singlet conformal operators present in this theory. Next, in section 2 
we derive supersymmetric Ward identities and resulting constrains on anomalous dimensions of our 
conformal operators. Then, in section 3 we proceed with the determination of universal non-forward 
anomalous dimension for Wilson twist-2 operators in N = 4 SYM. And, finally, section 4 contains 
our conclusion. 



Different supersymmetric Yang-Mills theories in four dimension can be constructed from higher 
dimensional N — 1 supersymmetric Yang-Mills theory through dimensional reduction [33] • To obtain 
the lagrangian of N = 4 SYM in 4-dimensions we start with N — 1 SYM in 10-dimensions 



Here V ^ are 10-dimensional Dirac gamma-matrices, G a v is gauge field strength, is covariant 
derivative and A a is Majorana-Weyl spinor. The supersymmetry transformations are 



1.1 Lagrangian of N = 4 SUSY YM 



£ = - G a p,G atu ' + % -\ a T^\ a . 



(1) 




a 



5 Q X 



(2) 



where = ifT^, r„]. After dimensional reduction to four dimensions we get the following la- 
grangian of iV = 4 supersymmetric Yang-Mills theory [33] 



£ = -± G ^G a ^+ l -\ a ^D»\ a + ±(D,A a r ) 2 + ±(D,B?) 2 

- 9 -f abc X a (a r A b r + l5 (3 r B b r ) X c (3) 

2 

- !L (f abc A b r A c t ) 2 + (f abc B b B') 2 + 2 (f abc A b r B') 2 , 



where A a denotes vector of 4 Majorana spinors and 




G' 



a 



a 



d,\ a + gf abc Al\ c , 




see last reference in [6'2\ for review 
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cr, 



^ = ( - 1 o ) • A=(f £ 

1,1 1 1 

7 [7m> 7f] > °>i = 7 K, a t ] , cr^ = - [#., ft] , n rt = - {a r , /3 f } 



In 4-dimensions we have gauge field (gluon), four Majorana fermions (gauginos), three scalars (A) 
and three pseudoscalars (B). Our notation is somewhat different form usually employed SU(4)- 
covariant form of N = 4 SYM, where we have gauge field, four left Weyl fermions and 6 real scalars. 
Certainly, both these forms of iV = 4 lagrangian are completely equivalent and we have chosen this 
one only because it seems to be more convenient in doing actual loop calculations. The 4-dimensional 
supersymmetry transformation (J2J) now become (here £ stands for a vector of 4 Majorana spinors) 

5AI = 

5A^. = £a r X a , 

5B a r = h,(3 r \\ (4) 
5X a = [a tlu G a ^+i lll D^{a r A a r + 1 ^ T B a r ) 

+ gf abc (a' rt A b r A c t + a Tt B b r B c t + 2 lh K rt A b r Bf)] £ , 

6X a = U-tvG^+zV^r^-Ts/W) 

- gf abc [a' rt A b r A\ + a rt B b r B c t - 2 l5 K rt A b B^)] . 

Note, if we choose A3 = A4 = 0, Ai, A2 non-vanishing and A\ = A2 = B\ = B2 = 0, A3 = A, B3 = B 
non-vanishing, then using the explicit representation of a and f3 matrices lagrangian (j3J) becomes 
lagrangian for N = 2 SUSY theory. 

Moreover, N = 4 SYM has internal 577(4) i?-symmetry group, which is symmetry of lagrangian 
under the following transformations of fields 

8A« = 0, 

5A a r = A' rt A« + A rt B?, 

5B a t = A rt B a r -A rt A a r7 (5) 

SX a = --(a' rt A' rt + a r tA rt + 2^ 5 k H A r t)X a , 

5X a = ^X a {a' rt A' r t + °rtA r t-2 l5 k rt A rt ), 

where parameters of these transformations, given by real antisymmetric matrices A^, A rt and A rt . 

Later in this paper we will be interested in scale transformation properties of Wilson operators 
of lowest twist or, what is the same, in those, which have maximal Lorentz spin. The component 
with maximal Lorentz spin has a symmetric traceless Lorentz structure and the simplest way to 
project it out is to use a convolution with the product of light-like vectors n M (n 2 = 0). This way we 
effectively project our theory on the light-cone. For a general four-vector X^ we introduce light-cone 
coordinates as follows: X ± = 1/ V / 2(X° ± X 3 ) and AT*% = X+Y~ + X~Y + - X l Y\ where % = 1, 2. 
For spinors the appropriate projectors are A ± = ^7 ± 7 T A, so that A = A + + A~. At this step it is 
also convenient to fix gauge condition. We employ light-cone gauge defined by setting A a+ = 0. 
Then taking £ = £ + (so that 5 A + = 0) the restricted light-cone supersymmetry transformations 
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(J3J) become (7^ = 7 M — %i7- — n*"f+, where and project "plus" and "minus" components 
respectively) 



S Q A a± = .^+ 7 ± A a- ; 5 Q A a = £ + \ a -, 5 Q B a = e+ 75 A a ~, 

5 Q A a - = - 7 - 7 ^ + ^V + ^"5 + (A a + 755 a )e + , (6) 

In light-cone gauge the restricted supersymmetry transformations form the off-shell supersymmetry 
algebra. They are linear and form a closed algebra on the projected +- components of fields, defined 
as components having maximal spin. 

1.2 Conformal twist-2 operators in N = 4 SUSY YM 

Now, let us introduce the local singlet (with respect to internal SU(4)-symmetry group) conformal 
Wilson twist-2 operators appearing in this models for unpolarized and polarized cases (T51 El 1221 EH] 

iG _ rial, t a a \l-ir< 5 / 2 f^ >+ \ „± na± 



Of, = Gli(id + y-^\(^)ei u Gli, (8) 



0% = G^(id + y-^^)g^Gti, (7) 

o), = Vcf (^) \v , (9) 

o), = \x%^d + ) l 1+1 ,cf 2 (^]\\\ (10) 



where = A + iB is a complex scalar field, T> — d — d, d = d + d, g^ u = g^ u — n^n* v — n^n*, 



& u = md^C'ilC^)^, (ID 



9aii9(iv n *p n v an d G^(z) are Gegenbauer polynomials 
r w z) _ (-1)"2" r(n + i/) r(n + 2z.) /2 d» ^n+,-1/21 ^ 

In what follows we will restrict ourselves to the analysis of singlet conformal operators, just intro- 
duced, and tensor operator to be introduced in next section. However, in general, N = 4 SYM 
has much richer content of twist-2 conformal operators, sitting in different representations of S77(4)- 
group. Besides bosonic operators in other irreducible representations of 5'[/(4)-group, we can write 
down fermionic (by quantum numbers) operators formed by scalar-gluino and gluon-gluino fields. 
While these fermionic operators were already present in theories with less supersymmetry like N = 1 
SYM (HI 122] and Wess-Zumino model {36], here we encounter new type of operator - vector operator 
formed by scalar and gluon fields. Under restricted supersymmetry transformations the conformal 
twist-2 operators form a closed operator supermultiplet. A general procedure of constructing super- 
multiplets of conformal operators is known for a long time jH] . Recently, full set of twist-2 conformal 
operators together with their transformations under restricted susy transformations was derived for 
the case of iV = 4 SYM in Ref. [37] and we refer interested reader to that paper. 
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As was already mentioned in introduction, we will be interested in the renormalization properties 
of these operators. It should be noted, that in the singlet case there is mixing between bosonic 
operators formed by gluon, gluino and scalar fields (|7|)- (fTT|) . Also, in the non-forward kinematics, in 
contrast to forward case, the operators (f7|)- (jll|) will mix under renormalization not only with each 
other, but also with the total derivatives of themselves. 

2 Supersymmetric Ward Identity in TV = 4 SUSY YM 

To begin with, let us introduce the following multiplicatively renormalized combinations of conformal 
operators for unpolarized (Eqs. (J7J), and (JTTJ) ) and polarized (Eqs. (JB1) and (fTUjO cases 

s h = QO h + J 2°i,i + J -^i^ t^ (13) 

Vjj = 6dl + 3 -d^, (14) 
Sh = 6^-^- (j + 1 ^ + 2) ^, (15) 

V l = 6 ^-4^' (16) 
Sit = 6^-^^+ (j + 2) 4 (j + 3) ^, (17) 

where the coefficients in front of operators can be found in a way, similar to Ej. Here, we 
would like to note, that these linear combinations of conformal operators, which are also the com- 
ponents of N = 4 operator supermultiplet, renormalize multiplicatively only in LO of perturbation 
theory. Beyond LO in dimensional reduction, which preserves supersymmetry to rather high order 
in perturbation theory, used in this paper they get additional rotations due to the breakdown of 
super conformal symmetry. So, in general, it is only the constrains on anomalous dimensions of our 
conformal operators following from supersymmetric Ward identities, that remain valid to all orders 
of perturbation theory. 

To derive the supersymmetric Ward identities, relating anomalous dimensions of singlet conformal 
operators to the anomalous dimension of some supersymmetry primary operator, which renormalized 
multiplicatively, we need to know the action of four restricted supersymmetry transformations on 
these singlet operators. Applying four restricted supersymmetry transformations to our initial singlet 



operators we get 2 : (5 Q = 5% '<§ '6°) 

S Q Sj tl = \(l-(-l) j %lvTte2liTtiW^ 2>l , (18) 

s Q sh = ^-(-ly^V^Vei wjj, (19) 

S Q Sli = ^(1 - (-1)^)^7-^367^7^1 W^, (20) 
S Q Vl = ^(l + (-l) J )67,V667,V6^, (21) 
ggj = ^(1 + (-1)^)^7-^7-6 , (22) 



2 Supersymmetric transformations for all operators in N = 4 SYM can be found in Appendix of Ref. 37 
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where 



IT 



3-1 



2d 



Gi 



(23) 
(24) 



r 



1 ( „Wr,va 



g u ± + g v f — g 1 ^ g p °) ■ Note, that here we kept only terms containing at the end 
operator we are interested in. 

We see one general property, which is inherent to this set of transformation: for bosonic operators 
with a certain parity (S % or V 1 ) the index j of final operator Wj^ changes on two units after each 
step. 

Using restricted supersymmetry transformation above, it is easy to find [22J, that the renormal- 
ized supersymmetric Ward identity in the regularization scheme, preserving supersymmetry, has the 
following form (Sjj denotes vector of operators Sh, Sh and S^) 



([S fl ]S Q X) = -(S Q [Sji]X) - (i[S j i]{6 Q S)X) and (S Q [S 3l }X) = finite. 



(25) 



where we used the fact, that renormalized action in supersymmetric regularization is invariant with 
respect to supersymmetry transformations (i[Sji](5 Q S)X) = 0. As we already noted, beyond LO 
our operators 1)130 . f|15f) and (|17|) mix under renormalization, so that the renormalized operators are 
defined as (square brackets correspond to renormalized quantities) 

\ / S 1 \ 

(26) 

kl 

and as a consequence the renormalization group equation for these operators is given by 



' s 1 ' 






s 2 






s 3 




k=0 \ 


















Z7 







)i 









d 



rfln/i 



" s 1 ' 




s 2 




s 3 





^5 1^5 13 7 5 \ 




' s 1 ' 


%yS 22^S 23^S 


1 


s 2 


XyS 32^S 33 7 <S J 


jk 


s 3 



(27) 



ki 



Now, from supersymmetric Ward identity (|25|) we get (a k = |(1 - (-l) fe ) and Z jk = for k > j) 



EE 



k=0 k'=0 

1/e poles in (j2~%|) cancel, provided 



lly 12y 12 y 
AS 

21y 227 237 

31y 32y 33y 



Ok 



{Zw}k-2,k> 
{■Zyy }k,k' 
{%W }k+2,k' 



finite . 



:28i 



E 

fc=0 

3 



117W 



fc-2J 



E 

fc=0 

3 



E{ 21 4 ] }.^[>v fc -d+E 

fc=0 J fc=0 



12 7 [1] 



3 k 



13y[l] 



k=0 
J 



3 k 



<T k [W, 



k+2,l\ 



22 zf) ^[W fe ,]+V{ 23 4 1] ) a k [W, 

) jk V J jk 



k+2,l\ 



3^2 \ Z W 

k=0 

3 ^2 { z w 



32 7 [1] 



fe=0 
J 



3 k 



33 7 [1] 



fc=0 
3 



3-%k 

[mi 



ki\ 



3,k 



k=0 

7 



3 k 



k+2,l\ 



<*3 



E * 



^[i] 



k=0 



j+2,k 



[W kl \. 



Taking into account linear independence of operators [Wjw] we finally get the following relations on 
anomalous dimensions of conformal operators 



11 S , 12 S j_ 13 S _ W 

72n+5,2m+l 72n+5,2m-l ~t~ 72n+5,2m-3 — 72n+3,2m-l) 



and 



21 S 22 $ 23 «5 

72ra+3,2m+l + 72n+3,2m-l + 72n+3,2m-3 

31 -/ + 32 -Y 5 + 33 -Y S 

/2n+l,2m+l ^ /2n+l,2m-l ' /2n+l,2m-3 



12_ 5 13 <S „ 

72n+l,2n+l i" 72n+l,2n-l = U ' 

22 ^ 5 + - ^ 

/2n+l,2n+l ' /2n+l,2n-l ~~ /2n+l,2n+l J 

72n+l,2n+l ^ /2n+l,2n-l — /2n+3,2n+l 5 



/2n+3,2m-l! 
/2n+3,2m-l5 



13 S 

72n+l,2n+l : 

/2n+l,2n+l 



m < n + 2 , 
m < n + 1 , 
m < n 



w 



72n+l,2n+l ~~ 72n+3,2n+3 ■ 



(29) 
(30) 
(31) 

(32) 
(33) 
(34) 



Now, let us turn to polarized case. All the steps one need to perform here are in one to one 
correspondence with unpolarized case. First we write down the supersymmetric Ward identity ("Pj,; 
denotes vector of operators Pj l (JH|) and Ph (HHJ)) 



{\Vji}5 Q X) = -{8 Q [Pji]X) - (i[V j i](5 Q P)X) and (S Q [V ] i]X) = finite. 



(35) 



The operators (fT4*|) and (fT6*j) mix under renormalization and thus we define renormalized operators 
as (square brackets correspond to renormalized quantities) 



V 1 
V 2 



jl 



E 

k=0 



llv 12y 

Zj-p Zj- 

21y 227 

Zi-p Zj- 



z; 1 
z- 1 



V 1 
V 2 



(36) 



A7 



The renormalization group equation for these operators is given by 



d 



V 1 
V 2 



E 



iXyP nf 

1\J> 22JP 



3 l k=0 

From supersymmetric Ward identity (|3*3)l we get ((Xfc = |(1 + 

j k 





' V 1 ' 


\k 


V 2 _ 



(37) 



■\) k ) and Z jk = for k > j) 



EE 

k=0 k'=0 

1/e poles in (jHH|) cancel, provided 



lly 12y 
Z/-p Zj 



21y 22y 
Zi'p Zj 



Ok 



E 

fc=0 

j 



u 7 [i] 

Zjj, 



E 

fc=0 



12 7 [1] 
Zjj, 



i Z W }k+l,k> 



[W, 



jk 



OkiW; 



k+l,l\ 



k'l\ 



3 

^E 



finite . 



(38) 



7W 

7v 



A;=0 

i 



j-l,fc 



A7 



fc=0 



fc=0 



fc=0 



Taking into account linear independence of operators [VV&i] we finally get the following relations 

m < n + 1 , (39) 
m < n (40) 



/2n+2,2m ^ /2n+2,2m-2 ~~ /2n,2m-2 ) 



21 P _|_ 22 V _ W 

/2rc,2m + 72n,2m-2 ~~ 72n,2m-2 J 
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and 



1 V - 

72n,2n — u 



22 >JP _ W 
I2n,2n ~ I2n,2n 



(41) 



Note, that both operators Wj^ and Wj^ are components of the same operator supermultiplet and 
thus have the same single anomalous dimension (with shifted argument). 

To obtain from these results the corresponding relations for anomalous dimensions of conformal 
operators considered earlier, one just need to use the following transition formulae among anomalous 
dimensions of conformal operators (JT3j) -(|17| ) and anomalous dimensions of operators (J7 | )-([li p : 



2\^P 22^P 




99^ 9>^Zj 




k 
2 

fc+3 
2 



-1 



(42) 



I 1^5 12^S 13 7 5 \ 

2y> 2^5 23 7 5 
\ 3^5 32^5 33^5 I 



( 



3(3+1) 
4 

(j+l)(j+2) 
12 



\ „ i+3 (j+2)(j+3) 

v 6 -— — i — 



/ 



/ 99- 
\ 



7 


9^ 




\ 


/ 6 


fc 
2 


k(k+l) 
4 


\ 


7 








6 


1 

4 


(fc+l)(fc+2) 
12 




7 




# 7 


/ 


\6 


fe+3 
2 


(fc+2)(fc+3) 
4 


/ 



(43) 



These relations allow us to operate both with anomalous dimensions of conformal operators (J7j)- ([TT 
or, equivalently, with anomalous dimensions of multiplicatively renormalized operators ()13j) - (jl7j) . 



3 Universal non-forward AD in N = 4 SUSY YM. 

Now, let us proceed with the determination of universal non-forward anomalous dimension 7™ = 
7 J^. The simplest way to do it is to use the formalism of broken conformal Ward identities (CWI) 
Ref. [T%1 119j. The basic idea of this method lies in the relation between scale and special conformal 
anomalies of conformal operators. It was found that non-diagonal part 7 ND of the complete anomalous 
dimensions matrix 

7*. = 7^,* + t}?, 7 = ^V> + (^) 2 7« + ... (44) 

arises entirely due to the violation of special conformal symmetry. Moreover, within this framework 
we have a relation between scale anomalous dimension matrix in n-th order of perturbation theory 
and matrices of scale and special conformal anomalies in (n— l)-order of perturbation theory. As a 
consequence the calculation of leading nondiagonal part of nonforward anomalous dimensions matrix, 
which is nonzero starting from two-loop order, could be reduced to calculation of more simple one- 
loop diagrams. 

Up to the moment there are results for nondiagonal parts of nonforward anomalous dimension 
matrices of conformal operators calculated in QED [20], QCD [2*T] . N — 1 supersymmetric Yang- 
Mills theory [22] and in supersymmetric Wess-Zumino model jHHl, both for polarized and unpolarized 
cases. The results in N = 1 SYM could be easily derived from already known QCD result via simple 
identification of colour Casimir operators: Ca = Cp = 2NfTp = N c . In N = 4 SYM the situation 
is slightly more complicated. Besides correspondingly adjusted map of colour Casimir operators: 
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Ca = Cf = \NfTp = N c , where we accounted for four Majorana fermions compared to one in 
N — 1 SYM, we need to consider scalar fields <p, which are absent in iV = 1 SYM. The presence 
of scalar fields results in introduction of one additional conformal operator (JT7j) and increases the 
size of anomalous dimension matrix in unpolarized case. So, in general, in unpolarized case we need 
to compute 5 new conformal anomalies. However, as we already noted, to find leading nondiagonal 
part of nonforward anomalous dimensions matrix all we need to know is one-loop matrix of special 
conformal anomalies (matrix of scale anomalies is diagonal at one-loop order and is known for a long 
time). A simple diagram analysis shows, that the only nontrivial contribution to special conformal 
anomaly matrix, not related to one-loop scale anomaly matrix and not already computed in QCD, 
resides in ^7 and ^7. 

On the contrary, in polarized case we have only two conformal operators: (JHJ) and (jlOjl. which 
are similar to polarized operators in iV = 1 SYM and QCD. In this case we do not need to per- 
form any new calculations and we can use polarized case to extract universal non-diagonal part 
of non-forward anomalous dimension matrix from already known results for non-forward polarized 
anomalous dimension matrix in QCD from Refs. fHl 1201 Ej and results for leading order forward 
polarized anomalous dimension matrix in N = 4 SYM from Refs. PUJ E]. As a result we obtain 
the following final expressions for leading non-diagonal parts of non-forward polarized anomalous 
dimension matrix in N = 4 SYM (j > k): 




(45) 



A ND(1) 




(46) 



|As~TO S! , _AA \g~V- 

+ Ij 9jk 9jk Ik 




(47) 



(48) 



where d 



■jk — 



Kj, k )/ a U, k ), 9jk 



Wjk/a(j, k) and 



a(j, k) 
b(j, k) 



2(i 
"[I 



-k)(j + k + 3), 
+ (_l)i-*](2A; + 3) 



if j > k and even ; otherwise . 



10 



The quantities Wj t k were computed in Refs. [T[JJ EE] and are given by 



'/// 



2[l + (-l)'-*]0 i _ 2 , fc (3 + 2A;) 



x 2Af A + (A jk - S x (j + 1)) 



u*o< 



+ 2(j-k)(j + k + 3) 



(*)< 



xx 



2[l + (-l)'- fc ]0 j _ 2jfc (3 + 2A;) 



x [2A jk + (A jk - S x (j + 1)) 



-2[l + (-iy-%_ 2 , fe (3 + 2A;) 







(j-k)(j + k + 3) 
{k + l)(k + 2) 
l{j-k)(j + k + 3) 



6 (A; + l)(fc + 2) 



where 9j k is equal to 1, if j > k and otherwise and 



A, 



Si 



j + k + 2 



Si 



j-k-2 



+ 2S 1 (j-k-l)-S 1 (j + l). 



(49) 

(50) 

(51) 
(52) 

(53) 



Here S k {n) denotes harmonic sum 



m=l 



and (a)^ = stands for Pochhammer symbol. The leading order forward anomalous dimensions 

of conformal operators, entering formula above are given by [TO! |B| 

1 



m=l 



-i)' 



m=l 



5i(m) (54) 



4 



(0) 



' 3 



u + m+2) 

48 

~(j + l)(j + 2) ' 
j(j+3) 
3(j + 1)(j + 2) ' 

45x0' + 1) + 



+ 5i0" + i) 



(7 + 1)0' + 2) ' 



(55) 
(56) 
(57) 
(58) 



Collecting all contributions together and using the relations (|4*2*|) and (jHSj) (or (|4*U|) ) we find the 
following expression for universal non-forward anomalous dimension in N = 4 SYM 



uni ND(1) 



w ^ (l + (-l)^)(2A;+l)(g 1 (j)-g 1 (fc)) 

7j-2,fc-2 — 4 ~ : ~ X 



Si 



j - k 



2S x (j-k)-S 



(k-l) 4 (j-k)(j + k + l) 
j + k 



U - 1) 4 + (fc - 1) 4 



+2(i-l) 4 5 1 (i)+2(fc-l) 4 5 1 (m. 



(59) 
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Now, combining Eqs. an d together we get 



iy + _ 2 t/ - 2 \ v 

I2n+2,2m ' l2n+2,2m-2 /2n,2m /2n,2m-2 



0. 



(60) 



Taking m = n — 1 it is easy to check using Eq. (|42j). that our nonforward anomalous dimensions satisfy 
this relation. Moreover we can check relation between the NLO non-forward anomalous dimension 
and NLO forward anomalous dimension from Ref. |7j 3 . Namely from Eq. ()60|) we have 



!y + 127 

/2n+2,2n ' /2n+2,2n-2 



2^ 
/2n,2n 



/2n,2n ^ /2n,2n-2 



/2n,2n-2 
/2n-2,2n-2 



o, 

0. 



for m 
for m 



n, 

n + 1. 



(61) 
(62) 



Indeed, taking explicit expressions for NLO forward anomalous dimensions from Ref. [7J (to adapt 
these results to normalization used in present paper one needs to multiply the results of Ref. [7j by 
(— §) , shift the value of momentum j by unity j = j + 1 and then make the following substitutions: 



h 6' 



22 V 



-4(S 3 (j) + SlaO')) + 85i(j)(5 2 (j) + S- 2 (j)) 
245x0-) 8(2j + 3) 



i(j' + l)(j' + 2) (./ • D 2 (./ • 2)'- 



(63) 
(64) 



where Si(j) stands for harmonic sum defined in Eq. ([54)1 . one can easily verify, that above relations 
hold true. Note, that Eq. (JoTj) contains a contribution from forward non-diagonal anomalous dimen- 
sion 21 72n,2n- The appearance of nondiagonal part of forward anomalous dimension matrix (}4*2*|) is 
related to the breaking of superconformal symmetry, what is explicitly demonstrated by Eq. Ijfilj) . 

Finally, having in mind possible future tests of AdS/CFT correspondence, we would like to 
consider one particular limit of our result Eq. (|59p. As this result is meaningful only for j > k we 
parameterize the limit of large j and k as follows: 



J 



nm 



k 



m . 



(65) 



Now, using S\(j) = *f?(j + 1) — ^(1) ~ hi(j) for large j, it is easy to find the following asymptotic 
behaviour of Eq. (|59j) (for even j — k) 



uni ND(1) 
In m,m 



16 ln(n) 



m rr 



n In 



n 



n- 



+ In 



rr 



(66) 



4 Conclusion 



In the present paper we found a closed analytical expression for NLO universal non-forward anoma- 
lous dimension of Wilson twist- 2 operators in N = 4 supersymmetric Yang-Mills theory. To derive 



3 In Ref. the results of calculation were written in DRED scheme, however the coupling constant was kept in 
MS scheme. To obtain a result fully expressed in DRED quantities, one needs to perform the following coupling 
constant redefinition for coupling constant ckdred = «ms + -^c/(127r) (ajyfg) 2 . It means, that we need to add to NLO 
anomalous dimension from Ref. |7| the corresponding LO expression from Refs. 10, 6 , that is in present paper we are 



using 7; 



„(!) 



Ref. [71 



1/3 7 



(0) 
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this result it was sufficient to know constrains on anomalous dimensions of conformal operators, 
following from supersymmetric Ward identities, together with already known results for anomalous 
dimensions of Wilson twist-2 operators in QCD. Moreover, supersymmetric Ward identities allowed 
us to find the relation of non-forward anomalous dimension determined here with already known 
forward anomalous dimension of Wilson twist-2 operators in iV = 4 SYM theory. This relation may 
serve as an additional check of correctness of the results obtained both for forward and non-forward 
anomalous dimensions of conformal operators in this model up to next-to-leading order in pertur- 
bation theory. Finally, we considered the limit of our result for non-forward anomalous dimension, 
when momenta j and k turn to infinity, while their ratio is being fixed. In this limit non-forward 
anomalous dimension scales as 1/j, where j is Lorentz spin. 
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